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ABSTRACT 

Recent observations reveal that a cool disk may survive in the innermost 
stable circular orbit (ISCO) for some black hole X-ray binaries in the canonical 
low/hard state. The spectrum is characterized by a power law with a photon 
index F ~ 1.5 — 2.1 in the range of 2-10 keV and a weak disk component with 
temperature of ~ 0.2 keV. In this work, We revisit the formation of such a cool, 
optically thick, geometrically thin disk in the most inner region of black hole 
X-ray binaries at the low/hard state within the context of disk accretion fed by 
condensation of hot corona. By taking into account the cooling process associated 
with both Compton and conductive processes in a corona, and the irradiation of 
the hot corona to the disk, we calculate the structure of the corona. For viscosity 
parameter a = 0.2, it's found that the inner disk can exist for accretion rate 
ranging from M ~ 0.006 — O.OSMEdd? over which the electron temperatures of 
the corona are in the range of 1 — 5 x 10^ K producing the hard X-ray emission. 
We calculate the emergent spectra of the inner disk and corona for different mass 
accretion rates. The effect of viscosity parameter a and albedo a (a is defined as 
the energy ratio of reflected radiation from the surface of the thin disk to incident 
radiation upon it from the corona) to the emergent spectra are also presented. 
Our model is used to explain the recent observations of GX 339-4 and Cyg X-1, in 
which the thin disk may exist at ISCO region in the low/hard state at luminosity 
around a few percent of LEdd- It's found that the observed maximal effective 
temperature of the thermal component and the hard X-ray photon index F can 
be matched well by our model. 

Subject headings: Accretion, accretion disks — Black hole physics — X-rays: 
individual (GX 339-4, Cyg X-1) — X-rays: stars 
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1. Introduction 

Black hole X-ray binaries (BHXRBs) are binary star systems that are luminous in the 
X-ray part of the spectrum. These X-ray emissions are generally thought to be caused by 
one of the component stars being a black hole accreting matter from the companion. The 
property of these systems have been reviewed by several authors (Remillard & McClintock 
2006; Done et al. 2007; Gilfanov et al. 2010). It is well-known that various spectral states 
have been exhibited in BHXRBs. In particular, two basic X-ray spectral states are presented 
with a high/soft spectral state occurring at high luminosities and a low/hard spectral state 
occurring at low luminosities. Generally, at the high/soft state, with higher accretion rate, 
the accretion is dominantly via a standard thin disk extending to the innermost stable 
circular orbit (ISCO) (Pringle & Rees 1972; Shakura & Sunyaev 1973; Mitsuda et al. 1984; 
Frank et el. 2002); The ISCO is at 3-Rs (where -Rs = 2GM/c^, G is gravitational constant, 
c is light speed, and M is the central black hole mass) for non-rotating black holes and 
0.615-Rs for rotating black hole with limiting spin rate a^, = 0.9982, where a^, is specific 
angular momentum (Thorne 1974); With decline of the accretion rate to some transition 
rate, BHXRBs enter the low/hard spectral state, in which the standard thin disk is replaced 
by a hot, optically thin, geometrically thick advection dominated accretion flows (ADAF) in 
the most inner region around the black hole, i.e., the thin disk truncates at some radius off 
the ISCO (Rees et al. 1982; Narayan & Yi 1994, 1995a, b; Abramowicz et al. 1995; Narayan 
2005, 2008; Kato et al. 2008; Esin et al. 1997; Kawabata & Mineshige 2010; Qiao et al. 
2010; Zhang et al. 2010). 

However, the picture of truncated outer thin disk -|- inner ADAF model in the low/hard 
state is challenged by the recent observations of some BHXRBs. Reis et al. (2010) investigate 
a sample composed of eight BHXRBs in the low/hard state observed by XMM-Newton and 
Suzaku. Although the hard X-ray continuum is characterized by a power law with a photon 
index in the range of F ~ 1.5 — 2.1, a thermal component with a color temperature consistent 
with L cx is detected in all eight sources, meanwhile, broad iron Kq, fluorescence line is 
also observed in half of the sample. Both the L oc relation and the broad iron line 
profile suggest that a cool disk extends to ISCO. The observed thermal component is around 
0.2 keV. Observations to GX 339-4 are made by Swift and RXTE at 2007 May when the 
transition to the low/hard state occurs (Kalemci et al. 2007). Fits to the observed broad iron 
line profile with relativistic reflection model at luminosity 0.8%LEdd suggests that a cool 
disk resides in the region very close to ISCO. A thermal component with a disk temperature 
~ 0.165 keV at luminosity 0.8%LEdd is detected. The X-ray spectra are roughly characterized 
by a power law with a photon index F = 1.63;*^q o3 (Tomsick et al. 2008). 

A disk accretion model maintained by the recondensation of hot corona has been pro- 
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posed to understand the presence of the cool disk in the most inner region around the black 
hole in the low/hard state by Liu et al. (2006), Meyer et al. (2007). In this model, when the 
accretion rate decreases just below the transition rate between high/soft state and low/hard 
state, as a consequence of efficient evaporation, the thin disk truncates at a distance where 
the maximum evaporation rate occurs (Meyer et al. 2000; Liu et al. 1999, 2002). The 
remnant disk from the truncation radius inwards can exist steadily fed by the condensation 
of the hot corona/ADAF rather than be swallowed by the black hole within the viscous time 
scale (Liu et al. 2007; Taam et al. 2008; Liu et al. 2011). The size of the inner remnant 
disk is governed by the accretion rate. The geometry of the accretion flows are allocated 
as an inner disk and a much cooler outer disk, which are separated by an ADAF. With 
decrease of the accretion rate, the inner disk shrinks and eventually vanishes completely at 
a certain accretion rate. The dynamical interaction between the inner disk and corona is 
studied by Liu et al. (2006), Meyer et al. (2007) in which only the vertical conductive cool- 
ing of the corona is considered throughout the corona. Furthermore, the Compton cooling 
of the corona by the soft photos from the underlying thin disk is taken into account. The 
cooling of the corona is classified as Compton-dominated case and conduction-dominated 
case at different radius. For Compton-dominated case, conductive cooling is neglected, and 
for conduction-dominated case, Compton cooling is neglected (Liu et al. 2007). However, for 
a certain mass accretion rate, Compton cooling and conductive cooling are comparable. In 
this case, both the Compton cooling and conductive cooling should be considered together. 
The cooling of the corona is modified by Taam et al. (2008), in which the effect of conductive 
cooling to the corona is added to Compton-dominated case, the effect of Compton cooling 
to the corona is added to conduction-dominated case. 

In this paper, in order to get a radially continuous, self-consistent disk-corona model, we 
employ the result of Taam et al. (2008) and consider the irradiation of the hot corona to the 
disk (Liu et al. 2011). We calculate the structure of the corona and the emergent spectra 
of the inner disk and corona with different parameters of the model (In our calculation, 
contribution of the ADAF gas in the gap between the inner disk and outer disk, and the 
outer disk and corona to the spectra are not included). For the emission of inner disk, both 
the accretion luminosity contributed from disk accretion and the irradiation of the corona 
are included. Our model is used to explain the recent observations of black hole X-ray 
binaries GX 339-4 and Cyg X-1 in the low/hard spectral state. In section 2, the disk-corona 
model is briefly described. In Section 3, we show the numerical results of the structure and 
the emergent spectra of the inner disk-corona system. In Section 4, we compare the model 
predictions with observations. Our conclusions are presented in Section 5. 
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2. The model 

The disk-corona model adopted here, is based on the study of Taam et al. (2008) 
and Liu et al. (2011). It is assumed that an ADAF-like hot corona (described by the 
self-similar solution of Narayan & Yi 1995, with a = 0.2, /3 = 0.8) lies above a thin disk. 
The corona is heated by the viscous release of gravitational energy of accreted gas and 
cooled by vertical conduction and inverse Compton scattering of soft photons emitted by the 
underlying disk. For the Compton cooling of the corona, the upward hard photons escape 
from the corona directly; the downward hard photons are partially reflected and partially 
absorbed by the underlying optically thick disk. The absorbed photons are reprocessed in 
the optically thick disk and remitted as a thermal emission providing the soft photons for 
the Compton scattering in the corona. This procedure is iterated until a stable disk-corona 
system forms. In the vertical transition layer between the disk and corona, an equilibrium 
is established between the conductive flux from the upper corona, bremsstrahlung radiation, 
and vertical enthalpy flux. For a given distance from the black hole, a fraction of the disk gas 
is heated and evaporated to the corona when the conduction flux is too large to be radiated 
away. On the other hand, a certain amount of coronal gas is cooled down, condensing to 
the disk if the bremsstrahlung radiation is more efficient than the conduction. At accretion 
rates around a few percent of Eddington value, gas evaporates from the disk to the corona, 
the disk vanishes at around a few hundred Schwarzschild radii, and the coronal gas partially 
condenses back to the disk in the innermost region. 

The cooling of the corona is classified as conduction-dominated case and Compton- 
dominated case at different radius (Liu et al. 2007). Here, we follow the work of Taam et al. 
(2008) in which the conduction-dominated case is modified by adding the Compton cooling, 
and the Compton-dominated case is modified by adding the conductive cooling. Meanwhile, 
in our calculation the irradiation of the corona to disk is added for getting a self-consistent 
disk-corona system (Liu et al. 2011). Throughout the paper, we scale some quantities as, 
m is the central black hole mass scaled with solar mass Mq. m is the accretion rate scaled 
with Eddington accretion rate MEdd, ^Edd = 1-39 x lO^^mg s~^. r is the distance from the 
black hole scaled with Schwarzschild radius Rs, Rs = ^GM/c^ = 2.95 x lO^mcm. For the 
sake of clarity, we list the basic results of the disk-corona model in the following. 

2.1. Corona Dominated by Conductive Cooling 

In the case of conduction dominated cooling the condensation rate is given as (Taam et 
al. 2008) 
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where A = ^p^j^^.^ is the ratio of Compton coohng rate to conductive coohng rate, which weighs 
the modification of neglected coohng process compared to the dominated coohng process. 



^cff,max. where T^^, 



The ratio can be reexpressed as A = 1.4052 x lO^mr 1 - 
is the maximum effective temperature of the accretion disk (Taam et aL 2008). r\ and Tq 
represent the size of the inner disk without and with modification in conduction-dominated 
case respectively, ri and Tq are expressed as, 



ri = 0.815a"2^/^m^/^ 
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(3) 



The electron temperature distribution of the corona in radial direction is. 
Tee = 2.01 X 10^0a-2/5m2/5r~2/5(l + xy^^K. 



(4) 



The Bremsstrahlung luminosity from the transition layer is 
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and the Compton luminosity from the corona is 
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Note that the emission from the outer pure ADAF is not included here, it would cause 
the coronal luminosity deviating from the true value when the inner disk is very small, which 
occurs at accretion rates less than 0.01. 
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2.2. Corona Dominated by Compton Cooling 

In the case of Compton dominated cooling the condensation rate is given as, 
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where 
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and the condensation radius Tq meets the following formula, 
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The electron temperature distribution of the corona in radial direction is, 

T,^ = 3.025 X 109a-2/5m-2/5rh2/5^i/5 
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The expression of Bremsstrahlung luminosity from the transition layer is same as equa- 
tion (IS]), and the the Compton luminosity from the corona is, 

^ = 0.392«-^/5m3/5m7/5fZk=y'^' 
LEdd V0.3 keV/ 

X r^\i + \y^/^x~^'/'''ii-x-'/^)dx. (11) 



For given black hole mass m, accretion rate rh, viscosity parameter a, and a presumed 
value of Tcff max, quantities describing the inner disk and corona are determined by eqs. ([T]) 
to in the case of conduction dominated cooling or by ([7]) to f|TT]) in the case of Compton 
dominated cooling. Iterative calculations are carried out until a self-consistent Tcg^max is 
obtained. 
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2.3. Iterative Calculations with Irradiation 



The surface effective temperature of the accretion disk with irradiation from the corona 
is expressed as, 



end 



(12) 



where Fend refers to the flux which originates from disk accretion fed by condensation in per 
unit area, i.e.. 
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-Firr refers to the the illumination flux from the corona to the disk surface per unit area (Liu 
et al. 2011), 
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where the ADAF-like corona is assumed above the disk as a point source at a height Hs, and 
the covering factor of the point source to a disk ring at distance R is given as / = ^^^^2^^2)3/2 ; 
a is albedo. L^au is the total intrinsic luminosity of the corona and the transition layer above 
and below the disk. 

Combining eqs. (fT2l) - (fni) the effective temperature of the accretion disk can be re- 
expressed as. 
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Mcndc2 

Tp'jj j^g^^ refers to the maximum effective temperature from disk accretion which is reached at 
'"tmax = (49/12). The expression of T^f^^^^ is given as (Liu et al. 2007), 
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At lower accretion rate m, it's assumed that the corona is dominated by conductive 
coohng. An effective temperature Teg max is presumed to calculate the condensation rate 
from eqs. ([T]) and ([2]), the luminosity of the transition layer and the corona from eqs. 
and with which a new effective temperature Teg^max is derived from eqs. f lTB]) and ffT7|) 
by assuming a value of albedo a. An iteration is made till the presumed temperature is 
consistent with the derived value. From the derived effective temperature, the Compton- 
dominated region is determined the following equation (Liu et al. 2007), which is. 



If eq. f|T8|) has no solution, it means the corona is dominated by conductive cooling through- 
out the corona, we find a self-consistent solution of the disk-corona system. Otherwise, the 
Compton-dominated region is determined by eq. fITH]) . We recalculate the condensation rate 
by combing eqs. ([T]), ([2]) and ([7]), ([H]), and the luminosity of transition layer and corona 
from eqs. ([5]), ([6]) and ( ITTjl till the presumed temperature is consistent with the derived 
value, we find a solution of disk-corona system. With increase of the mass accretion rate, 
Compton-dominated region extends inward and outward till Compton cooling dominates 
throughout the corona at some accretion rate. An effective temperature Tefj max is assumed 
to calculate the condensation rate from eqs. (E]), (|H]), and the luminosity from eqs. (E]) and 
f lTTjl . with which a new effective temperature Teg^max is calculated from eqs. f|T6|) and f|T7|) . 
Iterations are made till the presumed temperature is consistent with the derived value, we 
find a solution of the disk-corona system. 



With determination of the disk and corona features, we are able to calculate the spec- 
trum of the disk and corona. The contribution to the emergent spectra are composed of two 
components: (1) Multi black-body emission of the underlying thin disk which is partially 
scattered by the electrons in the hot corona (Compton emission) (2) Bremsstrahlung emis- 
sion of the transition layer between the disk and corona. We use Monte Carlo simulation 
to calculate the Comptom spectra from the hot corona. The method of the Monte Carlo 
simulation is based on Pozdniakov, Sobol' & Sunyaev (1977). We assume the electrons in 
the corona have a Maxwellian distribution, and the Compton scattering is the main cooling 
process. Since the electron-scattering optical depth in the corona is less than one, we intro- 
duce the weight w as described by Pozdniakov et al. (1977) in order to efficiently calculate 
the effects of multiple scattering. We first set wq = 1 for a given soft photon, then calculate 
the escape probability Pq of passing through the slab. The quantity of woPq are the trans- 
mitted portion and is recorded to calculate the penetrated spectrum or reprocessed photons 




(18) 



2.4. Calculation of the Spectra 
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according to the escape direction of the photon. The remaining weight wi = wo{l — Po) is the 
portion that undergoes at least one scattering. If we write the escape probabihty after the 
n-th scattering as the quantity WnPn is the transmitted portion of photons after the n-th 
scattering, and is recorded as upward or downward transmitted spectrum. The remaining 
portion Wn{l — Pn) undergoes the (n + l)-th scattering. This calculation is continued until 
the weight w becomes sufficiently small. The whole process is simulated by the Monte Carlo 
method. 

Because the temperature of the transition layer is obviously higher than the underlying 
cool disk, we expect the inverse Compton scattering of the photons from the transition layer 
in the hot corona is not important Compared with the soft photons from the underlying 
cool disk. So, we ignore the inverse Compton scattering of the bremsstrahlung photons from 
the transition layer in the hot corona. The bremsstrahlung spectral flux from the transition 
layer per unit area can be expressed as, 

pbrcm _ ^ „■ 7 

= /^e-'^^/'^^^p' ergs s~^ cm-^Hz'^ (19) 

where jp refers to the volume emissivity coefficient which is, ju = 6.8 x 10~^^ Z^neniT^^^^"^ 
^(z/,Tcpi)e-''''/'=^-pi, and / = 3.4 x IQ-^^Z^rieniT^J^^Z^pi. Z^pi is the thickness of the tran- 
sition layer, and Tcpi = 1.98 x 10^a~^/^m^/^K is the electron temperature of the transition 
layer (Liu et al. 2007). In eq. (19), the coefficient 1/2 means that only the upward photons 
from the transition layer can directly escape from the disk-corona system. The frequency- 
integrated flux from the transition layer is. 
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Solving eq. (120|) . we can get. 



J, _ ^ brem'" - ^2^^ 

the expression of Fbrcm is given as (Liu et al. 2007), 

Fbrem = ^ X 6.391 X 10^^a-'^^'^m~^m^/\-^/^ ergs s"^ cm'^, (22) 
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where the integral J^gg^iy, Tcp])e~^dx is a frequency average of the velocity averaged Gaunt 
factor, which is in the range 1.1 to 1.5. Choosing a value of 1.2 will give an accuracy to 
within about 12% (Rybicki & Lightman 1979). 

By combing eq. (fT9|l and eq. (pT]) . we get the local emergent bremsstrahlung spectra 
from the transition layer, which is given as, 

F^-- = f I ^ ^(^,T.pOe-W^^-p- ergs s-^ cm-^Hz-\ (23) 

By integrating the local Compton spectra from the corona and the bremsstrahlung 
spectra from the transition layer along radial direction, we get the total emergent spectra of 
the inner disk-corona system. 



3. Numerical Results 

Given a black holes mass m, the mass accretion rate m, the viscosity parameter a and 
the albedo a, the structure of the corona in radial direction is solved self-consistently as 
described in Section [2j Throughout the calculation, m = 10 is adopted. 

For viscosity parameter a = 0.2, and albedo a = 0.15, it's found that the inner disk 
can survive for accretion rate in the range of m ~ 0.006 — 0.03, which are consistent with 
the results of Liu et al (2011), where the modification to the cooling of the corona are 
not considered. This is because, at the upper limit of the accretion rate, the corona is 
dominated by the Compton cooling, the modification of conductive cooling to the corona 
can be neglected, and at the lower limit of the accretion rate, the corona is dominated by the 
Conductive cooling, the modification of Compton cooling to the corona can be neglected. 
The size of the accretion disk as functions of mass accretion rate is plotted in the left panel of 
Fig. [T]with red line. The size of the inner disk is defined by a critical radius, where neither the 
evaporation nor the condensation precess occurs. Outside the critical radius, the matter is 
evaporated to the corona from the disk (if the disk is not evaporated completely) and inside 
the critical radius, the corona matter condenses onto the disk. An increase of the mass 
accretion rate leads to more gases condensation onto the disk and the corresponding critical 
radius of the inner disk increases. The maximum temperature of the inner disk is presented 
in the right panel of Fig. [T]with red line. It's evident that the maximum temperature of the 
disk increases with the mass accretion rate. This is because, an increased value of the mass 
accretion rate results in both the increase of the condensation rate and the luminosity of the 
hot corona, the maximum temperature of the inner disk heated up by accreting condensed 
gas and the irradiation of the hot corona integrated throughout the inner disk increases. The 
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ratio of luminosity dissipated in the corona to the luminosity dissipated in the disk, Lc/La-, 
as functions of mass accretion rates m is plotted in Fig. [2] with red line. With increase of 
the accretion rate m, the relatively quick increase of the emission from the accretion disk to 
the emission from the corona leads to a decreased value of Lc/L^. The electron temperature 
and the optical depth for the Compton scattering of the hot electron in the vertical direction 
as functions of radius are plotted in Fig. [3] with red line. The solid line is for fn = 0.03 
(dotted line m = 0.02; dashed line: m = 0.01). We can see the electron temperatures are in 
the range of ~ 1 — 5 x 10^ K for different accretion rate 0.01 < m < 0.03. It's clear that the 
temperature of the corona decreases with increase of the mass accretion rate for a = 0.15. 
This is because, with increase of the mass accretion rate, the Compton cooling of the soft 
photons from the disk to the corona becomes more efficient, the temperature of the corona 
decreases. To clearly show the mass distribution of the inner disk-corona system in radial 
direction, we plot the mass accretion rate in the accretion disk and the mass accretion rate 
in the corona as functions of radius in Fig. HI It can be seen that, with increase of the mass 
accretion rate, more gases in the corona recondense back to the disk, the accretion rate in 
the disk increases. 

The emergent spectra of inner disk and corona with mass accretion rate are plotted in 
Fig. [5l The X-ray spectrum in the range of 2 — 10 keV is roughly characterized by a power 
law with a photon index F = 1.63 at m = 0.01 (solid green line). With increase of the mass 
accretion rate, the hard X-ray photon index is F = 1.75 at m = 0.02 (solid blue line) and 
F = 1.92 at m = 0.03 (solid red line). It's obviously that the X-ray spectra become softer 
with increase of the mass accretion rate. This can be understood as, the hard X-rays are 
dominated by Bremsstrahlung at low accretion rate, whereas by Compton radiation at high 
accretion rate. In Fig. [5l the dashed line is the contribution of the Compton scattering of 
the soft photons from the underlying cool disk by the electrons in the hot corona (Compton 
radiation), and the dotted line is the contribution from the bremsstrahlung of the transition 
layer. At lower accretion rate rh = 0.01, the cooling of the corona is dominated by the 
conductive cooling, fewer soft photons from the underling thin disk are scattered to the 
hard X-ray band, so the Compton emission to the hard X-ray band is weak, the hard X-ray 
emission is dominated by the bremsstrahlung of the transition layer. With increase of the 
mass accretion rate to rh = 0.02, the Compton emission of the corona to the X-ray band 
becomes dominant compared to the bremsstrahlung of the transition layer, as shown by the 
blue dashed line and dotted line in Fig. [51 With further increase of accretion rate, the 
Compton ^/-parameter decreases, which results in a large photon index and small luminosity 
ratio between the corona and disk, as shown by the red solid line in Fig. O and Fig. [21 

In order to study the effect of albedo a on the spectra, we plot the emergent spectra 
of the inner disk and corona with albedo in Fig. [6l The size and the maximum effective 
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temperature of the inner accretion disk with a can also be seen in Fig [H The electron 
temperature and the optical depth for the Compton scattering of the hot electron in the 
vertical direction as functions of radius are plotted in Fig. Ei The red line, blue line and 
green line are for albedo a = 0.15, a = 0.6 and a = 1 (without irradiation) respectively. 
For viscosity parameter a = 0.2 and mass accretion rate m = 0.02, we plot the emergent 
spectra for different albedo a in Fig. [61 For albedo a = 0.15, the hard X-ray photon index 
r = 1.75. With increase the value of albedo, the hard X-ray photon index are F = 1.63, 1.53 
for a = 0.6, 1 respectively. It can be seen clearly. X-ray spectrum becomes harder with 
increase of albedo. This is because for larger albedo, fewer X-ray photons of the irradiation 
from the corona to the disk are absorbed and reprocessed as thermal seed photons to be 
scattered in the corona, the temperature of the electron in the corona is higher, a harder 
X-ray spectrum is predicted. We plot the emergent spectra with albedo a for different mass 
accretion rate in Fig. [71 We can see, for both m = 0.01 and rh = 0.03, the hard X-ray 
spectra become harder with increase of a. From Fig. 6, it's also very clear that the positive 
dependence of the hard X-ray spectrum index on the mass accretion rate holds for different 
albedo a. 

From theory, the value of albedo a is uncertain, and has been investigated by several 
authors, e.g. George & Fabian (1991), White et al. (1988). The value of albedo a, which 
depends on how the incident photons from the corona interact with the underlying accretion 
disk, is defined as the energy in the reflective spectrum divided by the energy in the incident 
spectrum. The value of albedo is energy-dependent. White et al. (1988) study the Compton 
scattering of the incident photons by the electrons in the disk. In their calculations. X-rays is 
assumed to incident upon a semi-infinite, plane-parallel, zero temperature, purely scattering 
medium. They show that the monochromatic albedo a(xo) changes from 95% to 1% for 
the energy xq in the range of ~ 10"'^ to ~ 30 ( Xq is the energy of the incident photons 
scaled with rest energy of electron). If the mean energy of the photons are around 50 keV, 
~ 40% of the irradiation flux is absorbed by the disk through the Compton scattering of the 
electron in the disk (White et al. 1988). The K-shell absorption of the metals is negligible 
for the incident photons at ~ 50 keV. However, because the K-shell absorption of the metals 
is proportional to (mgC^/Ziz/)^/^, at lower energy, the K-shell absorption becomes important. 
For the distribution of N^, oc u~^'^, ~ 30% of the incident energy is absorbed (Taam et al. 
2008). Thus, ~ 70% of the total irradiation flux is absorbed by the Compton scattering and 
K-shell absorption, so the value of albedo a is around 0.3. However, it's highly uncertain for 
the calculation the value of albedo, e.g. the ionization state of the surface of the accretion 
disk. In this paper, we only take albedo as a parameter to fit the observations. 

To test the effect of the viscosity parameters a to the shape of the spectra, we plot the 
emergent spectra of the inner disk and corona with viscosity parameter a in Fig. [HI In our 
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calculation, albedo a = 0.15 and mass accretion rate m = 0.03 are adopted respectively. For 
a = 0.2, the hard X-ray photons is F = 1.92. With increase of a, the hard X-ray photon 
index are F = 1.67,1.57 for a = 0.25,0.3 respectively. It's found that the luminosity of 
the disk-corona system decreases systematically with the value of a, meanwhile, the X-ray 
spectra in the range of 2-10 keV becomes harder with the value of a. This is because, 
for fixed mass accretion rate, the surface density in the corona decreases with increase of 
a. With a lower density, the heat flux to the transition layer decreases, meanwhile the 
bremsstrahlung cooling rate is decreased even more. The energy balance between heating 
and cooling results in a decreased heating rate associated with a lower enthalpy flux. The net 
effect results in a decreased condensation rate. The energy conversion efficiency of the inner 
disk + corona/ADAF can be expressed as 77/0. 1 = {L / L-^dd) / {M / MEdd) ■ The efficiency 
of the thin disk is 0.1, however, for the hot corona/ADAF, due to the strong advection 
effect, the efficiency is much less than 0.1. With increase of a, due to the decrease of the 
condensation rate, less matters are accreted in the form of thin disk, the luminosity of the 
inner disk -|- corona/ADAF system decreases. Meanwhile, due to relative lack photons from 
the underlying thin disk to be scattered in the corona, the temperature of the corona is 
higher, so a harder X-ray spectrum is predicted. 

To show the effect of the different parameters (accretion rate rh, viscosity parameter a 
and albedo a) to the condensation features and the corresponding spectral features of the 
inner disk and corona clearly, we list the numerical results for different parameters in Table 
1. 



4. Comparison To Observations 
4.1. GX 339-4 

The distance to GX 339-4 is dynamically constrained at 6.7 < d < 9.3 kpc (Zdziarski 
et al. 2004), which is consistent with the lower limit of d > 6 kpc for kinematical distance 
based on analysis of high-resolution spectra of the Na D line (Hynes et al. 2004). The 
mass function of GX 339-4 is constrained to ~ 6Mq (Hynes et al. 2003; Mufioz-Darias et 
al. 2008). In this paper, a likely distance of d = 8 kpc and mass m = 5.8 are adopted to 
fit the recent observations (Tomsick et al. 2008). A thermal component with an inner disk 
temperature of 0.165 keV at a lower luminosity L ^ 0.008LEdd (1-100 keV) are detected by 
the observations of Swift and RXTE during 2007 June 10-14. The most constraints of the 
disk come from the reflection component, which requires that the disk can extend to r < 5 
(Tomsick et al. 2008). 
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With m = 5.8, the X-ray luminosity L = O.OOSLEdd and effective temperature of the 
disk Teff^max = 0.165 from observations. We assume viscosity parameter a = 0.3, and find 
that the observed X-ray luminosity and the effective temperature are matched by adopting 
rh = 0.0456 and albedo a = 0.782 respectively. The fitting results are listed in Table 2. We 
calculate the emergent spectra of the inner disk-corona system with the fitting parameters, 
and plot the spectrum in Fig. [9] (the solid line). It's found that the hard X-ray spectrum 
in the range of 1 — 100 keV is roughly characterized by a power law with a photon index 
Tmod = 1-50 from our model, which is roughly close to the observed one T = 1.631q;q3 
(Tomsick et al. 2008). 

With increase of the viscosity parameter to a = 0.4, we find the observed X-ray luminos- 
ity and the maximum effective temperature of the disk can be matched by taking accretion 
rate m = 0.0773 and albedo a = 0.617 respectively. In this case, the condensation rate is 
'^cnd = 5.98 X lO^'* which is smaller than mend = 1.4 x 10""^ for the case a = 0.3. So from 
eqs. (fT6l) and (fT7|) . a smaller value of albedo a is needed to match the observed effective 
temperature. We plot the corresponding emergent spectrum in Fig [9] (the dotted line). The 
hard X-ray spectrum from 1 — 100 keV is also roughly characterized by a power law with a 
photon index T^od = 1-4 from our model. The detailed fitting results are shown in Table 2. 

From the calculation above, it's found that, the hard X-ray photon index calculated by 
taking a = 0.3 can match the observed one better than that of taking a = 0.4. Here, we want 
to show that the observed X-ray luminosity and the maximum effective temperature of the 
disk can only be fitted in a very narrow range for the value of a. For instance, for a = 0.2, 
we can't fit the X-ray luminosity and the maximum effective temperature simultaneously, 
i.e., by taking m = 0.0251, the observed X-ray luminosity can be matched, however, even by 
taking the limited value of albedo a = 1, the theoretical value of Teff^max is 0.1825 keV which 
is bigger than the observed value 0.165 keV. 



4.2. CygX-1 

Cyg X-1 is a well-known black hole X-ray binaries and has been well studied by several 
authors. The distance to Cyg X-1 is very early estimated at 2 kpc (Murdin & Webster 1971; 
Reis et al. 2010), and confirmed by Massey et al. (1995) where a distance of 2.1 ±0.1 kpc is 
found. The mass of Cyg X-1 is found in the range of 7 — 15 Mq (Shaposhnikov & Titarchuk 
2009). In this paper, a distance d = 2 kpc and m = 10 are adopted (Reis et al. 2010). We 
compare our result with the recent observations of Cyg X-1 (Reis et al. 2010). They find that 
an X-ray luminosity Lx/^Edd = 5 x 10~'^(m/10) x {d/2 kpc)^ in the range of 0.5 — 10 keV and 
a hard power-law photon index F = 1.71 ± 0.01 (2-10 keV). Extrapolating the luminosity 
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to 100 keV, we get a X-ray luminosity Lx/^Edd = 1-47 x 10"2(m/10) x (d/2 kpc)^. By 
fitting the spectra with Diskbb+Power law model, they reveal a cool disk component with a 
temperature kT fa 0.194keV, meanwhile it's found that the inner boundary of the accretion 
disk is very close to ISCO, i.e., r^^ = 5.7+^'° (Reis et al. 2010). 

With m = 10, the X-ray luminosity = 1.47 x 10~^LEdd and effective temperature 
of the disk kT 0.194keV derived from observations, and assuming viscosity parameter 
a = 0.3, it's found that the observed X-ray luminosity and effective temperature are matched 
by taking m = 0.044 and albedo a = 0.398 respectively. We plot the emergent spectrum of 
the inner disk and corona with the fitting parameters in Fig. [10] (the solid line). It's found 
that the hard X-ray spectra in the range of 2 — 10 keV is characterized by a power law with a 
photos index Fmod = 1-60 from our model. The fitting results for bigger viscosity parameter 
a = 0.4 are also listed in Table. 2. 



5. Conclusion 

We study the formation of a cool, optically thick, geometrically thin disk in the most 
inner region of black hole X-ray binaries in the low/hard spectral state within the context 
of disk accretion fed by condensation of hot corona. By taking into account the cooling 
process associated with both Compton scattering and vertical conduction in a corona, and 
the irradiation of the hot corona to the disk, we obtain a self-consistent solution of the 
inner disk and corona. We calculate the emergent spectra of the inner disk and corona for 
different accretion rate, and examine the effect of different parameters on the spectra (e.g. 
viscosity parameters a and albedo a). Our model is used to explain the spectral features 
of BHXRBs GX 339-4 and Cyg X-1 in the low/hard spectral state, in which the thin disk 
resides in the region very closed to the ISCO. Our results are roughly in agreement with the 
observations characterized by a weak disk component with temperature of ~ 0.2 keV and a 
hard power-law X-ray spectra with F ~ 1.5 — 2.1 for the sources in the low/hard spectral 
state. 

We thank Prof. Ronald E. Taam for his very useful suggestions and comments. We 
thank the referee for his/her expert comments and suggestions to our paper. This work is sup- 
ported by the the National Natural Science Foundation of China (grants 11033007, 11173029) 
and by the National Basic Research Program of China-973 Program 2009CB824800. 
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Fig. 1. — The left panel: The size of the inner accretion disk as functions of mass accretion 
rate rh. In our calculation, a = 0.2 is adopted. The red line is for albedo a = 0.15 (blue 
line: a = 0.6; green line: a = 1). The right panel: The maximum temperature of the disk as 
functions of mass accretion rate rh. 
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Fig. 2. — The ratio of Lc/L^ as functions of mass accretion rate rh. In our calculation, 
a = 0.2 is adopted. The red line is for albedo a = 0.15 (blue line: a = 0.6; green line: 
a = 1). 
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Fig. 3. — The left panel: the electron temperature of the corona as a function of radius. In 
our calculation, a = 0.2 is adopted. The red lines represent the results calculated for albedo 
a = 0.15 (blue lines: a = 0.6; green line: a = 1). The different line types represent different 
accretion rates (dashed line: m = 0.01; dotted line rh = 0.02; solid line: rh = 0.03). The 
right panel: the Thomson optical depth for the Compton scattering of the hot corona in the 
vertical direction as a function of radius. 
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Fig. 4. — The mass accretion rate in the accretion disk and the mass accretion rate in the 
corona as functions of radius. The red line and black line represent the mass accretion rate 
in the disk and the mass accretion rate in the corona respectively. The solid line, dotted line 
and dashed line are for rh = 0.03, rh = 0.02 and rh = 0.01 respectively. In our calculation, 
a = 0.2 and albedo a = 0.15 are adopted. 
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Fig. 5. — The emergent spectra of inner disk and corona for different mass accretion rate. 
In our calculation, a = 0.2 and albedo a = 0.15 are adopted. The red line, blue line and 
green line are the spectra for m = 0.03, rh = 0.02 and rh = 0.01 respectively. The solid 
line is the total spectrum, the dashed line is the contribution from the Comptonization of 
the soft photons of the thin disk by the electrons in the hot corona, the dotted line is the 
contribution of the bremsstrahlung from the transition layer. 
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Fig. 6. — The emergent spectra of inner disk and corona for different albedo a. In our 
calculation, rh = 0.02 and a = 0.2 are adopted respectively. The solid red line is the total 
spectrum for a = 0.15, the red dashed line is the contribution from the Comptonization of 
the soft photons of the thin disk by the hot corona, the red dotted line is the contribution 
of the bremsstrahlung from the transition layer. The solid blue line and green line are for 
a = 0.6 and a = 1 respectively. 
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Fig. 7. — The emergent spectra of inner disk and corona for different accretion rate and 
albedo. In our calculation, a = 0.2 is adopted. The solid line, dotted line and dashed line 
are for rh = 0.03, rh = 0.02, and m = 0.01 respectively. The red line, blue line and green 
line are for albedo a = 0.15, a = 0.6 and a = 1 respectively. 
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Fig. 8. — The emergent spectra of inner disk and corona for different viscosity parameter 
a. In our calculation, m = 0.03 and a = 0.15 are adopted. The red solid line is the total 
spectrum for a = 0.2, the red dashed line is the contribution from the Comptonization of 
the soft photons of the thin disk by the hot corona, the red dotted line is the contribution 
of the bremsstrahlung from the transition layer. The solid blue line and green line are for 
a = 0.25 and a = 0.3 respectively. 
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Fig. 9. — The emergent spectra of GX 339-4 calculated by different fitting parameters of 
our model. The solid line is the emergent spectrum for m = 5.8, a = 0.3, m = 0.0456 and 
a = 0.782; The dotted line is the emergent spectrum for m = 5.8, a = 0.4, m = 0.0773 and 
a = 0.617. 
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Table 1: Condensation and spectral features of the inner disk and corona around the black 



hole of 10 Mr 
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0.3 
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-3 
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Note. — With black hole mass m=10, viscosity parameters a, albedo a, and the mass 
accretion rate rh, the size of the inner disk r^, condensation rate mend integrated from the 
condensation radius to 3-Rs! the luminosity dissipated in the corona Lx{cor)/-^Edd5 the ratio 
of luminosity dissipated in the corona to the luminosity dissipated in the disk Lc/L^, the 
maximum temperature of the inner disk Teff^max and the hard X-ray photon index Fmod in 
the range of 2-10 keV are listed. 



Table 2: Fitting Results for GX 339-4 and Cyg X-1 



GX 339-4: (m = 5.8^ L^/L^aa = 0.008^ Teff,max = 0.165 keV^ T = 1.631°;°^'^) 



Tjl a Td WT-cnd -^x(cor) /-^Edd L^/L^ Tmod 

0.3 0.0456 0.782 349 1.40 x lO^^ 0.8 x 10^^ ]^ 

0.4 0.0773 0.617 16.05 5.98 x 10"^ 0.8 x lO'^ 6.59 1.40 



Cyg X-1: (m = 10", Ljl^^d = 0.0147^ T,s,nr.. = 0.194 keV^ T = 1.71 ± 0.01'^) 
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rh 
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^d 
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-^x(cor) / -^Edd 


Lc/ Ld 


Tmod 
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0.044 


0.398 


55.0 
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2.57 


1.60 


O.4.... 
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0.267 


30.4 
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1.47 X 10"2 
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1.52 



Note. — With the black hole mass m, the X-ray luminosity ivx(cor) /-Z^Edd and the effec- 
tive temperatur Tcfr^max of the accretion disk derived from observation, and assuming a 
value of a, two quantities accretion rate rh and albedo a are adjusted simultaneously to 
match the observed X-ray luminosity and the effective temperatur of the accretion disk, rj, 
"^cnd, -^x(cor) /-^Edd and L^/ are the corresponding size of the inner disk, the integrated 
condensation rate from the condensation radius to 3-Rs, the ratio of luminosity dissipated 
in the corona to the luminosity dissipated in the disk respectively. Fmod is the hard X-ray 
photon index in the range of 2 — 10 keV calculated by the corresponding parameters. 

aThe central black hole mass 

hThe observed X-ray luminosity 
f cThe observed temperature of the accretion disk 
3iThe observed X-ray photon index in the range of 2-10 keV 
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Fig. 10. — The emergent spectra of Cyg X-1 calculated by different fitting parameters of 
our model. The solid line is the emergent spectrum for m = 10, a = 0.3, m = 0.044 and 
a = 0.398; The dotted line is the emergent spectrum for m = 10, a = 0.4, m = 0.064 and 
a = 0.267. 



